ON THE NUMBER OF ZEROS OF CERTAIN RATIONAL HARMONIC 

FUNCTIONS 



DMITRY KHAVINSON AND GENEVRA NEUMANN 

Abstract. Extending a result from |KS 03| . we show that the rational harmonic function 
r(z) — z, where r(z) is a rational function of degree n > 1, has no more than 5n — 5 complex 
zeros. Applications to gravitational lensing are discussed. In particular, this result settles 
a conjecture |Rh Oil concerning the maximum number of lensed images due to an n-point 
gravitational lens. 



1. Introduction 

A. Wilmshurst |Wil 98j showed that there is an upper bound on the number of zeros of a 
harmonic polynomial f(z) = p(z) — q(z), where p and q are analytic polynomials of differ- 
ent degree, answering the question of T. Sheil-Small |SS 92j . Let n = deg p > deg q = m. 
Wilmshurst showed that n 2 is a sharp upper bound when m = n — 1 and conjectured that the 
upper bound is actually m(m — 1) + 3n — 2. D. Khavinson and G. Swiatek KS 03 showed 
that Wilmshurst 's conjecture holds for the case n > l,m = 1 using methods from complex 
dynamics. When hearing of this result, P. Poggi-Corradini asked whether this approach can be 
extended to the case f(z) = p(z)/q(z) — z, where p and q are analytic polynomials. 

In this note, we apply the approach from |KS 03| to prove 

Theorem 1. Let r(z) = p(z)/q(z) be a rational function where p and q are relatively prime, 
analytic polynomials and such that n = deg r = max(deg p, deg q) > 1. Then 

#{z € C :Hz)=z} < 5n- 5 

We note that the zeros of r(z) — z are isolated, because each zero is also a fixed point of 

Q(z) = r(r(z)), an analytic rational function of degree n 2 . This also follows from a result of P. 
Davis | Da 74] (Chapter 14) concerning the Schwarz functions of analytic curves. (The Schwarz 
function is an analytic function S(z) that gives the equation of a curve in the form z = S(z), 
cf. |Da 74j .) A rational Schwarz function implies that the curve is a line or a circle, so the 
degree must be one. 

We also note that r{z) — z will not have a zero at oo. If oo were a zero, then r(z) = 
az + b + 0(1/ z). By the change of variable z = 1/w, we see that a/w — 1/w + b + 0(w) has 
a zero at w = 0. Restricting w to the real axis, we see that a = 1 and b = 0. We obtain a 
contradiction when we restrict w to the imaginary axis. 
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We now discuss an application of our result to gravitational microlensing. An n-point grav- 
itational lens can be modeled as follows. Suppose that we have n point masses (such as stars). 
Construct a plane through the center of mass of these point masses, such that the line of sight 
from the observer to the center of mass is orthogonal to this plane. This plane is called the 
lens plane (or deflector plane). Suppose that the lens plane is between the observer and a light 
source. (We are assuming that the distance between the point masses is small compared to 
the distance between the observer and the lens plane, as well as the distance between the lens 
plane and the light source.) The plane containing our light source which is parallel to the lens 
plane is called the source plane. Due to the deflection of light by the point masses, multiple 
images of the light source are formed. This phenomenon is known as gravitational microlens- 
ing. See jWa 98 and NB 96 for an introduction to gravitational lensing and |St 97j for an 
introduction to a complex formulation of lensing theory; also see |PLW Ol] , 

Gravitational microlensing can be modeled by a lens equation, which defines a mapping from 
the lens plane to the source plane. To set up a lens equation for our n-point gravitational lens, 
the point masses are projected onto the lens plane. The projection of the j-th point mass has 
a scaled mass of rrij and is located at a scaled coordinate of Zj in the lens plane, where rrij is 
a positive constant and Zj is a complex constant. Suppose that we have a light source located 
at a scaled coordinate of w in the source plane. Following |Wit 90j . this lens equation will be 
given by w = z + "f~z — sign(<j)X™ =1 rrij/ (~z — zj), where the normalized (external) shear 7 and 
the optical depth (or normalized surface density) a ^ are real constants. In this model, if z 
satisfies the lens equation, then our gravitational lens will map z to w; hence z corresponds to 
the position of a lensed image. The number of lensed images is the number of solutions of the 
lens equation. 

We can rewrite this lens equation in terms of the rational harmonic function f(z) = r(z) — z 
by letting r(z) = w — jz + sign(cr)S" =1 rrij/(z — Zj). We thus see that the zeros of f(z) are 
solutions of the lens equation for a light source at position w. H. Witt |Wit 90| showed for 
n > 1 that the maximum number of observed images is at most n 2 + 1 when 7 = and (n + l) 2 
when 7 7^ 0. A. Petters |Pe 92j used Morse theory to obtain further estimates for both cases. 
S. Mao, A. Petters, and H. Witt jMPW 97] conjectured that the maximum number of images 
is linear in n for the case 7 = and a > 0. S. H. Rhie |Rh 01] later conjectured that for n > 1 
such a gravitational lens gives at most 5n — 5 images. In the 7 = case, deg r = n and Theorem 
1 settles this conjecture. Further, for the case 7 7^ 0, we see that deg r = n + 1, so Theorem 1 
gives an upper bound of 5(n + 1) — 5 = 5n lensed images. As a result, we have the following 
corollary: 

Corollary 1. An n-point gravitational lens modeled by w = z +72 — s^fi(c)S™ =1 mj/(z — ~z~j), 
where n > 1 can produce at most 5n — 5 images when the shear 7 = and at most 5n images 
when the shear 7 7^ 0. 
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3. Preliminaries 

We first recall some terminology and a fact from complex dynamics. Let r(z) = p(z)/q(z) be 
a rational function, where p and q are relatively prime, analytic polynomials. The degree of r, 
denoted by deg r, is given by max(deg p, deg q). z is called a critical point of r if r'(z) = 0. We 
will be interested in counting the fixed points of a rational function. A fixed point zq £ Coo is 
said to be attractive, repelling or neutral if |r'(zo)| < 1, |r'(zo)| > 1, or |r'(zo)| = 1 respectively. 
A neutral fixed point where the derivative is a root of unity is said to be rationally neutral. A 
fixed point zq is said to attract some point w G if the iterates of r at w converge to zq. We 
will be using the following result, whose proof can be found in |CG 93j . Chapter III, Theorems 
2.2 and 2.3: 

Fact 1. Let r be a rational function with deg r > 1. If zq is an attracting or rationally neutral 
fixed point, then z$ attracts some critical point of r. 

We will also use a version of the argument principle for harmonic functions. Let / be harmonic 
in an open set. We say that z is in the critical set of / if the Jacobian of / vanishes at z. In 
particular, in a sufficiently small disk containing z, there exist analytic functions h and g such 
that / = h + g in this disk. Then z is in the critical set of / if |/i'(z)| 2 — |</(2i)| 2 = 0. We 
say that z is a singular zero of / if z is in the critical set of / and f{z) = 0. / is said to be 
sense-preserving at z if > |</(.z)| and is said to be sense-reversing at z if |^'(^)| < |</(z)|. 

The order of a non-singular zero is positive if / is sense-preserving at the zero and negative if 
/ is sense-reversing. 

Let / be harmonic in a punctured neighborhood of zq. We will refer to zq as a pole of / 
if f(z) — > oo as z — > zq. Let C be an oriented closed curve that contains neither zeros nor 
poles of /. The notation Ac arg f{z) denotes the increment in the argument of f{z) along C. 
Following |ST 00| . the order of a pole of / is given by — t^Ac arg f(z), where C is a sufficiently 
small circle around the pole. We note that if / is sense-reversing on a sufficiently small circle 
around the pole, then the order of the pole will be negative. We will use the following version 
of the argument principle which is taken from |ST OOj : 

Fact 2. Let f be harmonic, except for a finite number of poles, in a simply connected domain 
D in the complex plane. Let C be a Jordan curve contained in D not passing through a pole 
or a zero, and let £1 be the open bounded region created by C. Suppose that f has no singular 
zeros in D and let N be the sum of the orders of the zeros of f in CI. Let M be the sum of the 
orders of the poles of f in $7. Then Ac arg f{z) = 2nN — 2nM . 

We note that a more general form of the argument principle can be found in |SS 02| . 

4. NON-REPELLING FIXED POINTS 

Proposition 1. Let p and q be relatively prime analytic polynomials. If r = p/q is a rational 
function of degree n > 1, then the set of points for which z = r(z) and \r'(z)\ < 1 has cardinality 
at most 2n — 2. 
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Proof. Let n+ denote the number of points zq satisfying the conditions of Proposition ^ 
Following the approach of D. Khavinson and G. Swiatek |KS 03j . we consider the function 

Q(z) = r(r(z)), which is a rational function of degree n 2 . As in KS_03j, we note that all fixed 
points of r{z) that are critical points for f(z) = r(z) — z are rationally neutral fixed points 
for Q(z), so Fact ^ will a Pply- Their outline (Lemmas 1 - 3 in KS 03 together with Fact ^) 
carry over word for word to the rational case if C is replaced by the Riemann sphere Coo. in 
particular, each point zq which satisfies the conditions of Proposition ^ attracts at least n + 1 
critical points of Q. 

Since deg Q = n 2 , by the Riemann-Hurwitz relation (see |Fo 81j . Section 17.14), Q has 
2n 2 — 2 critical points (counted with multiplicities) in Coo- n + will be largest when all of the 
critical points of Q are attracted to points zq satisfying the conditions of Proposition ^ Since 
different fixed points attract disjoint sets of critical points, we see that 2n 2 — 2 > (n + 1) 
and the claim follows. □ 

5. Proof of Theorem [T] 

As in |KS 03| . we say that the rational function r(z) is regular if no critical point of f{z) = 
r{z) — z is a zero of /. This will allow us to apply Fact |2] to count the zeros of /. 

Ansatz 1. If r is regular of degree n > 1, then f(z) = r(z) — z has at most 5n — 5 distinct 
finite zeros . 

Proof. Let r = p/q, where p and q are relatively prime, analytic polynomials. Let n p = deg p 
and n q = deg q. Then n = m&x{n p ,n q }. Let ra_ denote the number of sense-reversing zeros of 
/ and let n + be defined as in the proof of Proposition 1. Since r is regular, no zero of / lies 
on the critical set of /; this will allow us to apply Fact [2 In particular, n+ counts the number 
of sense-preserving zeros of /. Then / has n + + n_ zeros in the finite complex plane, counting 
multiplicity. Note that / is sense-reversing at each of the n q poles of r; hence, the order of each 
pole is negative. 

Consider the increment in the argument of / in a region bounded by a circle of large enough 
radius so that all of the finite zeros and poles of / are enclosed. We will apply Fact [21 and use 
our bounds on n+ from Proposition ^ Note that we are counting zeros with multiplicity. We 
consider two cases. 

(1) Assume that n = n q > n p > 0. Then the critical set of / is bounded. For z large, r(z) is 
at most 0(1). Thus / is sense-preserving on our circle with an argument change of 2tt. 
We also note that oo is not a zero of /. ByFact|^l 1 = (n + — n_) — (— n) < 2n— 2— n^+n. 
Hence n_ < 3n — 3, so / has at most 5n — 5 zeros in Coo and the claim follows. 

(2) Assume that n = n p > n q + 1. Since p and q are relatively prime, at least one of the 
two polynomials has a non-zero constant term. 

We first suppose that p has a non-zero constant term. Then /(0) / 0. Let z = 1/w 
and consider F(w) = l/r(l/w) —w = G(w)/H(w) — w. Then F satisfies the conditions 
of the previous case, replacing n q by deg H and n p by deg G. Thus, F has at most 
5n — 5 zeros. Since /(0) 7^ 0, / must also have at most 5n — 5 zeros in the finite complex 
plane. 

We now suppose that p does not have a non-zero constant term. Since p and q are 
relatively prime, the lowest order term in q must be a non-zero constant. Consider 
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f c (z) = (p( z ) + c)/q(z) — z = r c {l/w) — 1/w. For all c sufficiently small, we have that 
p{z) + c and q are relatively prime. In that case, / and f c have the same poles and 
f c approaches / uniformly as c —* on compact subsets of C which do not contain 
any of the poles of /. We then substitute z = 1/w to form F c (w) = l/r c {\/w) — w 
= G c (w)/H c (w) — w. By construction, deg G c = deg H c = n, so F c has at most 5n — 5 
zeros in the finite plane. Hence, this also holds for f c since / c (0) ^ 0. 

Now suppose that / has more than 5n — 5 zeros in C. By constructing a sufficiently 
small circle around each zero of /, we can guarantee that / is harmonic in each of the 
resulting closed disks. Moreover, we can make our circles sufficiently small so that each 
closed disk contains no other zeros of / and contains no critical points of / (recall that 
r is regular). Assume that f c does not vanish in such a disk surrounding a zero, so 
|/ c | > a > 0. Then, for c sufficiently small, arg / = arg (f c (1 + (/ — / c )// c )) has zero 
increment around the circle bounding this disk, by Fact this contradicts / having a 
zero in this disk. We see that f c must have the same number of zeros as / in that disk 
(counting multiplicity). Thus, f c must have more than 5n — 5 zeros in the finite plane, 
a contradiction. □ 

It remains to show that it is enough to consider regular rational functions. 

Lemma 1. // r(z) is a rational function of degree greater than 1, then the set of complex 
numbers c for which r(z) — c is regular is open and dense in C. 

Proof. As in |KS 03] . it is enough to show that the image of the critical set of f(z) = r(z) — z is 
nowhere dense in C. In contrast to the polynomial case, it is possible for the critical set to be 
unbounded. This difficulty however is easily resolved by restricting / to an increasing family 
of concentric balls whose radii run to oo. □ 

Lemma Q shows that the set of regular rational functions is dense in the topology of uniform 
convergence in the spherical metric. As at the end of the proof of the Ansatz, it is easily seen 
that every zero of the function f(z) = r(z) — z must be a limit point for zeros of r{z) — z — c 
where c — > 0. Thus, Lemman an d the Ansatz prove the theorem. 

6. Final remarks 

Our proof does not indicate whether the 5n — 5 bound is always sharp. We have, however, 
found an example where this bound is attained for the case n = 2; namely, f{z) = (~z 2 + z — 
7})/(z 2 — ^z + 1) — z. This function has 5n — 5 = 5 distinct finite zeros. This function has 
two poles: z = ^(3 — i y/7) and z = |(3 + i y/7); f is sense-reversing at these poles. There are 
three sense-reversing zeros: z = |, i(l + i vTE), and |(1 — i \/lT). As expected from Fact |2] 
for an overall index change of +1 on large circles, there are two sense-preserving zeros; namely, 
z = 1 — \pl and 1 + \/2, Figure ^ is a Mathematica plot of the critical set of this function and 
Figure |21 shows its image. We note that / is sense-preserving in the unbounded component in 
Figure sense-reversing in the larger of the two bounded components, and sense-preserving in 
the smaller of the two bounded components. We also note that / cannot be rewritten to model 
a 2-point gravitational lens. 

L. Geyer |Ge 03] has recently shown that the 3n — 2 bound on the number of zeros of 
f(z) = p{z) — 1 where deg p = n is sharp for all n > 1. D. Bshouty and A. Lyzzaik |BL 03] 
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Figure 1. Critical set for f(z) 
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FIGURE 2. Image of critical set for f(z) 
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have recently given an elementary proof for n = 4,5, 6, 8. Hence, a sharp bound on the number 
of zeros of f(z) = r{z) — z must be at least 3n — 2. 

Results from gravitational lensing help with the question of sharpness of the 5n — 5 bound 
in Theorem 1. We restrict our attention to the 7 = 0, a > case of the lens equation (c/. 
Introduction). Mao, Petters, and Witt MPW 97] have shown there are 3n + 1 images when 
the light source is at the origin and the n > 2 point masses (each of mass 1/n) are vertices of a 
regular polygon centered at the origin. Thus, in contrast to the case where r{z) is a polynomial, 
by moving the poles of r(z) into the finite plane, we see that a sharp bound on the number of 
zeros must be at least 3n + 1. Further, by a clever perturbation argument, Rhie Rh 03] has 
shown that the 5n — 5 bound is attained and is hence sharp for all n > 1. 

It is known that the number of lensed images is odd in the case of gravitational lensing by 
a regular gravitational lens ( |St 97j . |Bu 81j ). given that the position of the source is not a 
critical value of the lens equation. (A regular lens has a smooth mass distribution. A point w 
is a critical value of the lens equation if the Jacobian of the lens mapping defined by the lens 
equation vanishes at any of the points in Given that 5n — 5 is not always odd, it may 

seem surprising that 5n — 5 lensed images may be possible for an n-point gravitational lens with 
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7 = and a > 0. However, it was shown by Petters |Pe 92j using Morse theory, that this is not 
a problem. We summarize this result in the following corollary and present an alternate proof. 

Corollary 2. Suppose that we have an n-point gravitational lens modeled as above with 7 = 
and a light source which is not located at a critical value of the lens equation. Then the number 
of lensed images must be even when n is odd and odd when n is even. 

Proof. For this case of the lens equation, we may apply step (1) in the proof of the Ansatz. 
In particular, the number of images N will be n+ + where n + denotes the number of 
sense-preserving solutions and n_ the number of sense-reversing solutions of the lens equation. 
From step (1), we see that n+ = 1 + n_ — n, so N = 1 + 2n_ — n. □ 

This result has also been shown by other techniques in |Rh Olj (footnote 2, page 2), which 
extends the approach of W. Burke |Bu 81j for the case of a regular gravitational lens to that of 
an n-point gravitational lens. Note that the above proof is similar to the approach used by N. 
Straumann St 92j for the case of a regular gravitational lens. 

We also note that Theorem 1 can be applied to the case of a more general mass distribution 
than point masses. For a compactly supported mass distribution with the projected mass den- 
sity dtp in the lens plane, the lens equation transforms into w = z + yz — sign(a) J c dip(() / (z — 
C) (cf. NB 96j and |St 97j ). In particular, we have the following corollary: 

Corollary 3. Suppose that the projected mass density of a gravitational lens consists of n > 1 
radially symmetric, continuous, compactly supported densities in the lens plane. If the shear 
7 = 0, then the number of lensed images outside the support of the masses cannot exceed 5n — 5. 
If the shear 7 / 0, then the number of lensed images outside the support of the masses cannot 
exceed 5n. 

Proof. Consider the integral term in the lens equation for one of the mass densities. A cal- 
culation shows that this integral evaluated for any z outside of the support of this radially 
symmetric mass density equals b/{z — a), where o is the center of the mass density and b is a 
finite constant. Hence, for z outside of the support of the masses, the lens equation reduces 
to the lens equation for n point masses (each mass is reduced to a point mass at the sphere's 
center). Corollary 1 can then be applied to count the possible images, giving the upper bounds 
on images outside the support of the masses as claimed. □ 

Remark 1. If the mass distribution is composed entirely of luminous matter, then the only 
lensed images that will be visible are those that lie outside of the support of the masses. However, 
there are mass distributions that are not luminous (see [APOD 0T| and |AP()D 03] ). For such 
mass distributions, it is possible to have lensed images that lie inside the support of the masses. 

A question related to bounding the number of zeros of f(z) = r{z) — z would be to find a 
bound on the number of zeros of f{z) = R(z) — r(z), where both R and r are rational functions. 
It is not clear what condition would be needed for such a function to have a finite number of 
zeros, much less a bound on the number of distinct zeros (beyond the obvious Bezout Theorem's 
bound when we assume that there are a finite number of zeros). For example, let p(z) be a 
polynomial (of degree at least two) and consider f(z) = p(z) — l/p(z). The zeros of this 
function form a lemniscate in the complex plane. In other words, f(z) has an infinite number 
of zeros and these zeros are not isolated. Moreover, lmx s _- >QO f{z) = 00. If we compare this with 
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the result of Wilmshurst |Wil 98| for the case when / is entire: f(z) — > oo as z — > oo implies a 
finite number of zeros, we see that all of the zeros for an entire function will be isolated under 
Wilmshurst's condition for finite valence. Thus, / having finite poles makes the problem of 
discreteness of its zero set much more subtle. 
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